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Consider a nonlinear system, with memory, which has two input ports 
and one output port. It is assumed that the system can be represented by 
a double Volterra series. Two results for such a system are stated in Part I. 
The first is a general expression for the sinusoidal components of the 
output y(t) when the two inputs x u (t) and x v (t) are sums of sinusoidal 
terms. The second result is an expression for the power spectrum of y (t) 
when x u (t) is a stationary Gaussian process and x v (t) = P cos pt. Part 
II is concerned with using results from the theory of Volterra series for 
multi-input systems to calculate the third-order distortion in an idealized 
frequency converter. 

I. INTRODUCTION 

This paper deals with nonlinear, time-invariant systems with 
memory which (i) have more than one input port, and (ii) are driven 
by inputs which are essentially sums of sine waves. 

The paper consists of two parts. Part I is concerned with a system 
which has two inputs, x u (t) and x v (t), and one output y(t). Two results 
for single-input systems are generalized : (i) an expression is given for 
an arbitrary frequency component of y(t) when x u (t) and x v (t) are 
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finite sums of sine waves, and (it) an expression is given for the power 
spectrum of y (t) when x u (t) is a stationary, zero-mean, Gaussian noise, 
and x v (t) is a single sine wave. 

Although most of the discussion in Part I deals with systems having 
two input ports, many of the results can be formally generalized to 
systems with more than two inputs. 

Part II is devoted to an example which shows how results given in 
Ref. 1 for a one-input Volterra system can be used to examine systems 
consisting of a single two-terminal nonlinear element imbedded in a 
linear network containing sources. The transformation from a multi- 
input to a single-input system, is based upon Thevenin's theorem 
(see, for example, Anderson and Leon 2 ). The example treated here 
is a frequency converter using a nonlinear capacitor. Particular atten- 
tion is paid to computing the limiting form of the expression for the 
third-order distortion when the signal and pump amplitudes become 
small. 

The procedure we use in Part II is essentially a systemization of a 
procedure used by Gardiner and Ghobrial 3 to study the distortion 
performance of a varactor frequency converter. As they point out, 
their treatment differs from the linear time-varying analysis usually 
employed to study frequency converters. It is appropriate to mention 
here that a promising new general method of computing distortion in 
frequency converters has been developed by R. B. Swerdlow. 4 His 
method is based upon the use of Volterra series with time-varying 
kernels. 

Part I. Two Input Ports 

When analysis of the type used to study Volterra systems is applied 
to nonlinear circuits having two input ports and one output port, 
some of the simpler results for one-input circuits can be generalized 
in a straightforward way. Here we state two such generalizations. In 
the first, the two inputs are sums of sine waves. In the second, one 
input is stationary, zero-mean Gaussian noise, and the other input is a 
single sine wave. 

The derivations of the generalizations are not given here because 
they consist of rather straigntforward, although lengthy, applications 
of the procedures used in Ref. 1 to deal with the one-input case. 

II. DOUBLE VOLTERRA SERIES 

Let x u (t) and x v (t) be the two inputs, and let the output y(t) be 
given by the double Volterra series 
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1/(0 = f f ' ^rr, /* x **■ • ■ r **■ r *i- ■ ■ r 



(2» , 



•g m ;n(Ul, ■ • -, U m ; Vi, ■ • •, V n ) II X u (t - U r ) fl ».(< ~ ».), (1) 

r=l «=1 

where the prime on £' means that the term m = n = is omitted. 
The product II is understood to have the value 1 when the number of 
factors (m or n) is 0, and if n, say, is zero there are no ^-integrations. 
The kernel g m; „ is a symmetric function of u h • • • , w m and of Vi, • • • , w„. 
For the inputs that we shall consider, the (m + n)-fold Fourier 
transform 



Gm;n(ful, ' - ', fum ', fvl, ' ' ', fi-n) = I tfttl' • • I 



dV n 



■g m -.n(ui, ■••;•• •«„) exp [— j(unaui H h *>„a> vn )] (2) 

plays an important role. Here G 0; o = 0, co = 2wf, and (? m!B is a sym- 
metric function of f uh ■ ■ ■ , f um and of f vh ■■■, /„„. 

Much as in Ref. 1, the "harmonic input" method can be used to 
determine G m -„ from the system equations by setting 

m 
*■(<) = L exp (j'wurO, 

'V (3) 

3.(0 = £ exp (iw„0, 

where the co's are incommensurable, and solving for the coefficient of 
exp Zj(co u i + • • • + w um +«,!+•••+ ai rn )i] in the expansion of y(0- 
This coefficient is equal to G m ., n (f ul , • • •, f um ; f vl , ■ ■■, /„„). Note that 
if the system output y(t) is applied to the input of a linear transducer, 
the transducer output can also be expressed as a double Volterra series. 
The transducer output function corresponding to the transducer input 
function G m;n is 

FlJ(o>ui + • • • + a>„ n )](? m;n (/ ul , ••■;•• ./,»), 

where F(jta) is the transfer function of the transducer. 

If, say, n is zero and m > 0, G m;0 (/ U i, / M 2, ■ ■ • , fun,;) is equal to 
the coefficient of exp [j("«i + w„ 2 + • ■ • + u um )Q in the expansion of 
y(t) when x v (t) = and x u (t) is given by (3). 

There is a resemblance between the double Volterra series (1) for 
the two-port inputs x u (t), x,.(t) (Case A) and the single Volterra series 
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for the special input x(t) = x*(t) + x v (t) (Case B). For Case B, 
y(t) = Z n / dui--- / du k g k (ui, •'•,uh) 

f[ [X u {t - U T ) + X v (t - «,)]. (4) 



r=l 



When the product is expanded and the symmetry of |* is used, the 
product can be written as 

E ( J Ux u (t-u r ) n *.(<-«*+.). 

Setting n = fc - m and w m+ . = ». for 8 = 1, 2, • • • , n carries (4) into a 
form which goes into (1) when g m+n (ui, • • •, u m , Vi, • • -,v n ) is replaced 

by ffm;»(Wl, •",«»; 01, '•-, *>»)• 

The results stated below for Case A show a similar resemblance to 
the corresponding Case B stated in Ref. 1. For example, when x u (t) 
and x v (t) are the sums of sinusoidal terms, the expression for a particu- 
lar component in y(t) for Case A can be obtained from the correspond- 
ing expression for Case B by replacing G m+n in Case B by G m -, n and 
inserting a semicolon at the appropriate place in the string of argu- 
ments [as in eq. (9) below]. 

III. SINUSOIDAL INPUTS 

When the inputs x u (t) and x v (t) are sums of sinusoidal waves, an 
expression for any particular component in y(t) can be obtained by 
extending the analysis given in Section VI-B of Ref. 1. There, eqs. 
[1, (138), (139), (140)] [meaning eqs. (138), (139), and (140) of 
Ref. 1] show that if the input x(t) to a one-input system is given by 

X(t) = t, Pr COS U>rt, (5) 



r=l 



where the oj r 's are incommensurable, then the exp [j(iVicui + • • • 
4- N^M, N r ^ 0, component of y(t) is 

« m r (p r /2) Nr+2lr 1 
exp[i(Ar 1 . 1 +...+ ^ M )CE o ---E o n[ (iVr + Zr)! ^ ! J 

■G n £(fih 1+h , (-/i)*„ • • -, (/mW^ (-/m)*J> (6) 
where G = 0, (/„)* denotes the string of fc arguments /„, /„, ■••, /#, 
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and the subscript n on G has the value 

n = t (N r + 2l r ). (7) 

r— 1 

Here the notation of Ref. 1 has been changed to bring the statement 
of (6) in line with the notation used in the present paper. 

Methods of computing (6) when the G n 's are constants, i.e., are 
independent of frequency, have been considered by several writers 
(see Kroupa 5 and Sea and Vacroux 6 ). 

To state the generalization of (6) let 

M x 

«u(0 = Y,Pr COS 0)4, X v (t) = £ PrCOSWr*, (8) 

r-1 r-M+1 

where the u> r 's are incommensurable, oi r = 27r/ r , X = /x + ", and n and p 

are positive integers. Then the exp Zj(Nia>i -\ h N\<a\)t], N r ^ 0, 

component in y(t) is 

exp [j(iVi Wl +. -.+ 2W*] £".£ IT [ 7^X7^ 1 

1,-0 J X =0r = l L (^r + tr) ! if ! J 

Gm;nL(fl)lfi+lv ( — fl)h, (/OlTtfln •••, (/m)jV m +I mI ("~fn)l/, 

• (A+OaVh i+Im+i> *'"' (h)Nx+h, ( — f\)hl- (9) 

Here G 0; o = 0, and if Z or iV + Z are the corresponding arguments do 
not appear in G m;n . The values of m and n are 

m = t (Nr + 2J r ), n = £ (# r + %.)• (10) 

r-l r= M +l 

The semicolon in the subscript of G m;n differs in meaning from the 
semicolon used in [1, (139), (140)]. The notation (/,)* is the same as 
that in (6) and in [1, (169)]. The series (9) may either converge or 
diverge, depending on the P's and G's. 

Changing the signs of a>i and /i in (9) carries (9) into the expression 

for the exp [j'(- Nim + N 2 oa H h Nxco^Q component in y(t), 

etc. [see the discussion below eq. (5) in Ref. 1]. When some of the 
wr's are commensurable, some of the components in y(t) coalesce and 
can be treated by the method used in [1, (6), (7)]. 

To examine the case in which x*(t) contains a dc component, let /i 
and wi tend to in (8) and (9). Then P x is the dc component of x u (t) 

and the exp [j(N 2 co2 -\ \- Nxw^f] component of y(t) is the result 

of the coalescence (as /i— >0) of (i) the components exp [_j(Nim 
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+ AT 2 a> 2 + ••• + NwM for tfi = 0, 1, 2, •••, oo and (n) 

exp [j'(- Nim + # 2 oj 2 H + ATxo, x )C for Ni = 1, 2, ■ • •, oo. When 

(9) and (9) with -/i in place of /i are summed over the values of Ni, 
the double sum with respect to h and Ni can be reduced to a single 
sum by setting k = Ni + 2h and using the binomial theorem. The 
desired component, namely exp [j (iV 2 «2 +•••+ N\(a\)f\, in y(t) 

when x u (t) = Pi + P 2 cos w 2 H hP, cos «„* and s,(0 is given by 

(8) is found to be 

„ 00 =0 pk x r (P r /2) Nr+2lr I 

exp U(Nw + • • • + ffx«x)<] E o E o - • • E o ^ n [ (A r r r + z r )u r! J 
•(^•[(O)*, (AW* (-/Oii. •■•;■••, (AVx+ix. (-/x)uJ (11) 

Here n is given by (10), and m = fc + (#2 + 2Z 2 ) H h (#p + 2Z„) 

when fj. = 2 and m = k when p = 1. 

Equation (9) can be generalized to the case of three or more input 
ports in a straightforward way. 

IV. X u (t) GAUSSIAN AND X v (t) = P COS pt 

The case x v (t) = P cos p« and x u (t) = I(t), where /(*) is a station- 
ary, zero-mean, Gaussian noise having the two-sided power spectrum 
WiU), can be handled in much the same way as was the case x(t) 
= 7(0 + P cos pt discussed in Section VII-C of Ref. 1. 

The discrete sinusoidal components in y(t) are given by the ensemble 
average 

(y(t)) = t c n exp(jnpt), (12) 



7i=— X 



where 



- (P/2)"+mi q , n 



Sn,a,k(fl, ••-, fk) fp) = L, -JS; (/2H-*i»H-|»lLJl» Jll » /» 

,,=0 " I * 

-/:, /l, A, ••',/*; («»/p)»+l»h (- Sn/,),]. (13) 

Here 2tt/ p = p, s» = 1 for n = 0, s„ = - 1 f or n < 0, and as in (6), 
(s„f P ) a denotes a sequence of <t arguments, all equal to s„/ p . As ex- 
plained in connection with [1, (145)], Q,[WtUY\ denotes a y-fold 
integration with respect to fi, •■-, f, with limits ± « . The integrand 
is Wi(fi)- ■ -Wi(fl) times the function [in (13) the function is G] 
of fi, •••, fl represented by all of the terms lying to the right of 
QrlWitf')]. Qo[WiU')l denotes the identity operator. 
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The two-sided power spectrum of y(t) is 
W V U)= t \cn\*Hf - nf p ) 



, y Q£W,(m y , (f , , 



L 



(P/2) 



2<H-|n| 



bCff+lnDlo-!^"'''"^ 1 ' ■"» /*»/>) 



(14) 



Replacing G M . B by (? m+n in eq. (13) for 5 and substituting in (14) 
gives the expression [1, (175)] for W y (f) in the single input port case 
x(t) = I(t) + P cos pi. There is a corresponding similarity between 
the one input port formula [1, (16)] when the two-port expression 
(14) for W v (f) is written out. 

V. EXAMPLE — COMPUTATION OF G\-\ 

Consider the circuit shown in Fig. 1. The admittance H(f) is linear 
but the resistor R and capacitor C are nonlinear. The voltage across 
ft is 

<*I» + PH (15) 

and the capacitance of C depends upon the charge Q, the capacitance 
being a + bQ. The output of interest is the voltage y(t) across C : 

Q = (a + bQ)y, 
/« + /, = / = dQ/dt. 

The current /„ (t) is given by 



(16) 



Lit) = (" h{u)[x v {t - u) - y(t - u)]du, 

J —X 



(17) 



where h(u) is the Fourier transform of H(f). 
Elimination of /„ leads to the circuit equations 

X u = al u + 0/jj + y, 
dQ/dt = I u + / h(u)lx v (t - u) - y(t - u)]du, (18) 

Q = (a + 6<?)y. 

The G's corresponding to y can be obtained from (18) by the harmonic 
input method. In using this method it is convenient to work with the 
notation z k = exp (jW) where the w's are incommensurable. 

In order to get G liq (/ x ;) we set x u = z h x v = 0, y = dz x + higher 
harmonics, /„ = i x z x -\ , and Q = q& x -\ . The harmonic input 
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Fig. 1— Circuit with input voltages x„(0, ».(<), output voltage y(t), and nonlinear 
R and C. 

method states that <?i,o (/ij) is equal to c x . Substituting in (18) and 
equating coefficients of Zi gives 



1 = ah + C\, 
jwiqi = t'i - H(fi)ci, 
51 = ac\. 



(19) 



Solving for Ci gives Gi :0 (/ij). Similarly, starting with x u = and 
x v = z\ gives (j 0: i(;/i). The results are 






(20) 



where the subscript A means that / = /i is to be substituted in 
a = 27r/and in H(f). 

To get (?i;i(/i;/a) we set x u = Zi, x„ = z 2) and assume 

2/ = CiZi + C2Z2 + C12Z1Z2 + ■ • • 1 

iu = ii*i + *a*i + ii2ZiZ2 H , (21) 

Q = ?izi + (72Z2 + 9122^2 H • 

When (21) is substituted in the circuit equations (18), the coefficients 
of zi give the equations (19) and therefore Ci = OiMi'i)- Similarly, 
c 2 = Go-, 1(1/2). The coefficients of Z1Z2 give a set of equations which, 
upon solving for c 12 and using qi = ac h q 2 = ac 2) t 2 = — c 2 /a, i\ = ■ • • , 
give 



en 



_ f (P/<*)(H + j«a)/! - i(«i + w 2 )aa6 "1 _ (22) 

~ - ClC2 L (1 + off + j«aa) /l+/l J 

Replacing Cic 2 by Gi 5 o(/i;)Go : i(;/«) S ives the required expression for 
GvAfi-Ji) = C12. § 

To get OrMiJi ;) we start with x u = z x + 2*, ». = and agam make 
the substitutions (21) in the circuit equations (18). The coefficients of 
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2i2 2 give the same set of equations as before because x u and x v ap- 
pear only linearly in the circuit equations. We have C\ = Oi-o(fil), 
ii = C\[H(fi) + ju\a~], q\ = aci, and q 2 = ac 2 as before, but now 
c 2 = Gi ; o(/ 2 ;), it = c 2 [#(/ 2 ) + jw&l, and c u = 0%Mi, fa;). 
To sum up, we have 

GiMilU) = 2Gi ; o(/i;)Go ; i(;/ 2 ) 

X C ex P r ession in brackets (22)], (23) 

where Gi i0 and G 0; i are given by (20). Expressions for G 2:0 (/i, fa;) and 
GoM'Jh fa) are obtained by replacing the product Gi :0 G -,i in (23) by 
Gi-oGi-o and Go-.iQo-.i, respectively, and changing the bracket slightly. 

To get (?i ;2 (/i; fa, f^) we set x u = Si, x, = 2 2 + 3 and proceed 
along the lines used to get Gi-i, and so on. 

As an example of the use of (23), suppose that x u = Pi cos uit and 
x v = P 2 cos co 2 i. Then the exp Ci( w i ± »a)0 component in y is, 
from (9), 



exp [j(»i ± oj 2 )<] 



PiPj 



Gi;i(/i; ±/,) + 



(24) 



Similarly, the leading terms in the series for the components of fre- 
quency 2/i and 2/ 2 are given by <? 2i o(/i, /ij) and G t t i i(,7si /»), 
respectively. 

Part II. Analysis for a Simple Frequency Converter 

Here the circuit shown in Fig. 2 is used as an example to show how 
a multi-input system can sometimes be analyzed by the single-input 
formulas of Ref. 1. The output of interest is the voltage y(t) across the 
nonlinear capacitor C. TheVenin's theorem is used to replace the circuit 
of Fig. 2 by that of Fig. 3, and a recurrence relation is derived for the 
corresponding G„'s. The results are used to get an expression for the 
third-order distortion when Fig. 2 is regarded as the circuit for an 
up-converter. 
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Fig. 2 — Frequency converter with nonlinear capacitor C. 
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VI. REDUCTION OF A MULTIPLE-INPUT SYSTEM TO A SINGLE-INPUT 
SYSTEM 

The system considered here and in the following sections is shown 
in Fig. 2. The admittances #„(/), H v {f), H t (f) are linear and C is 
the nonlinear capacitor used in Section V. The charge on C is Q(t), the 
current i" = dQ/dt flows into C, the capacitance of C is a + bQ, and 
the voltage y(t) across C is related to Q(t) by 

Q = (a + bQ)y. (25) 

Po is a biasing dc voltage. 

The problem is to determine the components of y(t) when 

X u (t) = Pi cos o>ii + Pi cos uit, 
X„(t) = Pp COS Upt, 

and «i, w 2 , and w p are incommensurable. 

As far as y(t) is concerned, the analysis of Fig. 2 can be reduced to 
that of the simpler circuit shown in Fig. 3. To accomplish this we 
apply Thevenin's theorem to the portion of Fig. 2 lying to the left 
of the terminals of C. As far as the exp OW) components of y(t) and 
I(t) are concerned, this portion of the system can be replaced by an 
admittance #(/i) - H u (fi) + #*(/i) + Hi(fi) ^ series with the 
(open-circuit) voltage 






Similar consideration of the remaining components shows that I(t) 
and y(t) can be computed from the circuit of Fig. 3 in which 

H(f) = H u {f) + H v (f) + H l U), 
x{t) = poPo + P\P\ cos (oii + <pi) + piPi cos M + ^2) 

+ PpPp COS (oipt + <p p ), 

The equation for y(t), namely 

I _^ = r h{u )lx{t -u)- y(t - u)ldu, (28) 

at 1 — by J-x, 

where h(u) is the Fourier transform of H{f), can be obtained by 
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Fig. 3 — Theveiiin equivalent of Fig. 2. 



equating two expressions for /, the one on the left being / = dQ/dt 
in which Q = (a + bQ)y = ay/ (I — by). 

It is convenient to subtract out the dc components of x(t) and y(t) 
and apply the formulas of Ref. 1 to the portion y(t) of y(t) which tends 
to when x(t) — > 0, x(t) being the time-varying component of x(t). 
Therefore, in the system equation (28), we make the substitutions 



x(t) = x + x(t), 
y(t) = y + y(t), 



(29) 



where x is the dc value of x(t) and y is the (dc) value of y{t) when 
x{t) = x . From (27) x = p P , and substitution in (28) gives 
= H(0) (xq - y ). Assuming H(0) ^0 gives y = x = poPq. It 
should be noticed that the substitutions (29) are not strictly necessary 
because the dc component in x(t) could be handled (at the cost of 
more work) by the analogue of (11). 

Subtracting the result of substituting x and y Q in (28) from the 
result of substituting (29) in (28) and using 



aQ/o + y) _ 
1 - b(y + y) 1 - fa/c 



ay 



= a y 

" 1 -by 



(30) 



shows that the system to be analyzed by the single-input formulas of 
Ref. 1 is described by the equations 



d 
dt 



&y _ r 
i - bu J- 



h(u)\jt(t — u) — y(t — u)~]du, 



(3D 



X(t) = piPi COS (wit + ^i) + piP-i COS (a) 2 t + ip 2 ) 

+ p p P p cos (co p t + ifp), (32) 



where 



d = o/(l — byt 



6 = 6/(1 -by ). 



(33) 



Here x(t) and y(t) play the roles that x(t) and y(t) play in Ref. 1; 
and in the remainder of this paper the G n 's will refer to x(t) and fl(t). 
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VII. CALCULATION OF THE G»'s 

The functions Gi(/i), 0t(/i>/i)»* ' ' can be computed from (31) by 
the harmonic input method. A guide to the work is furnished by the 
resemblance of our problem to the one described by Fig. 3 of Ref. 1 
and eqs. [1, (42), (43), (106)]. Expanding the left side of (31) as 

i t ^ W W(0]' (34) 

at i=i 

carries (31) into the form of [1, (106)] except for the operator d/dt. A 
procedure similar to the one used to deal with [1, (106)] gives 

<?»(/i, h) = 'bG 1 (f 1 )G 1 (f 2 )K(f l + h), 
Gt(f h U, f») = ^G 1 (f l )G i (f2)G 1 (fz)K(f 1 + /, + /•) (35) 

■\_KUx + h) + K(fi + /.) + K{ji + h) + 3], 
and the recurrence relation 

<?»(/i, •■•,/») = \KUx+ — + toi&*Q9Uu- ••»/-)• (36) 

1=2 

The G^'s are the G n 's for the Volterra series for [y(t)y, and formulas 
for computing them are given in [1, (24) to (29)]. Equation (36) is a 
recurrence relation because Gg> can be expressed as the sum of products 
of (?i, G 2 , • • •, Gn-i- By starting with (?i(/i), the relation (36) can be 
used to compute G 2 , G 3 , • ■ ■ in succession. In the next section, (36) will 
be used to compute G 4 - 

VIII. COMPONENTS OF y(t) OF FREQUENCY /i + f p AND 2/i — f» + /p 

In this section we use (6) and the recurrence relation (36) to derive 
expressions for the exp [j(wi + «,)*] and exp [j(2a>i - a> 2 + «,)«] 
components of y(t) in Fig. 2 when (i) P h P 2 , P P are small, («) /i and 
U are nearly equal, and (Hi) H(/) is zero except for frequencies lying 
in narrow bands about the values 

0, U, f P , /«, (37) 

where /„ denotes the upper sideband frequency /i + /„. 

The component of frequency 2/i - f 2 + f P represents a typical 
third-order distortion product in an up-converter when /i and / 2 are 
signal frequencies, f p the pump frequency, and / u = /i + / p , /2 + f P 
the output frequencies. 
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Table I — Notation for Various Values of H(f) and K(f) 



Frequency, / 


H(f) 


K(f) 


0, /i - fi 

/i, ft, 2/i - h 

U, fi-f* + f P 

A + /,(-/«), 2/1 -/, + /, 

outside bands 


H 
H x 

k 






K x 

K p 

K u 

-2 



As mentioned in Section I, the procedure we use here can be regarded 
as a systemization of a method used by Gardiner and Ghobrial 3 to 
study the distortion performance of a varactor frequency converter. 
In our notation, the problem they solve is that of determining the 
exp [i(2a>i — W2 + w P )0 component of the charge Q(t) from the 
system equation 

V(t) = aQ + bQ 2 + r k(u)I(t - u)du. (38) 

J — X 

Here V{t) is the sum of three sine waves [just as £ is in (32)], / = dQ/dt, 
and k(u) is the Fourier transform of the linear impedance Z(f) in the 
TheVenin equivalent of the converter circuit. It can be shown from 
(38) that the G„'s corresponding to Q(t) can be determined by recur- 
rence from 

GtUi) = l/(a+ju>Z) /v 

- b 1 



Gnifl,- • -,fn) 



-[.- 



+ juZ\ 



/l+"-+/n 



Q*Uu- ■■,/.). 



(39) 



Now we return to our own problem. From the expression (32) for 
the input x (t) and the leading terms in the series (6) it follows that the 
ex P [i(2o>i — o>2 + o) p )Q and exp [i(wi + u P )Q components of y(t) 
are, respectively, 

exp {i[(2wi - U2 + a> P )t + 2<pi - <p 2 + <p p ]\ 

■ (p?p2PpP?P 2 P P /32)[(?4(/i, h, -h, /,)+•••]. (40) 

exp {j[(o)i + o)p)t + ip x + ip,}) (piPpP iPp/4) 

■[G,(U, /,) + (p^/8)G 4 (/ 1) /„ /„ -/,)+•••]• (41) 

Only one G 4 term appears in (41) because we shall assume that 
PiPi/ppPp and P2P2/ PpPp are small compared to one. 

The function G 2 is given by (35), and the remaining problem is to 
compute (? 4 from the formula obtained by setting n = 4 in the recur- 
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rence relation (36) : 

(?4(/i, /., /■, u) = iK(fi + u + n + /«) 

■f^W/n /■,/•./«). (42) 

1=2 

As explained in Ref. 1 in connection with eqs. [1, (24) to (29)], we have 
■fiQPVufafaM = (D(2)(3)(4), (43) 

!jG1 3 >(/i, /,, /a, /«) = (D(2)(34) + (1)(3)(24) + (D(4)(23) 

+ (2)(3)(14) + (2)(4)(13) + (3)(4)(12), (44) 

± OFUh /•, /a, /0 - (D(234) + (2) (134) + (3) (124) + (4) (123) 

+ (12) (34) + (13) (24) + (14) (23), (45) 

where we have written "(2)," for example, for Oi(fi), "(34)" for 
G 2 (f 3 , /«), "(234)" for G a (/2, /b, A), and so on. 

The next step is to compute the right-hand sides of (43), (44), (45) 
from the expressions (35) for G h G 2 , G 3 when only frequencies in the 
bands indicated by (37) are allowed to flow. To aid in this, we introduce 
the notation shown in Table I for the values of the function K(f) 
= - 2jcad/(H(f) +jwd) needed for the various GYs and G 3 's. We 
make the usual assumption that the admittance H (/) remains constant 
in each band. 

We consider first the G 4 (/i, /1. ~ h, f P ) in expression (40) for the 
exp [i(2wi - co 2 + ta p )f\ component of y(t). Equation (43) gives 

GFUu h, -U, fp) = 24G?(/x)G 1 (-/ 2 )G 1 (/ P ), (46) 

where, from (35), Gi(f) = [H/{H +i«d)]/. In eq. (44) for G[ Z) 
"(34)" now means G 2 (-/2, /,) and from the expression (35) for G 2 
and Table I we get 

(34) = &<?i(-/s)<?i(/p)ff(-/2 + fp) = 6Gi(-/ 2 )G 1 (/ p )(- 2). 

Similarly, " (24)" means G 2 (fi,f P ) and using the notation K(/i+/ P ) = K u 
gives 

(24) = bG l (f 1 )G 1 (f p )K u , 

and so on. Going through all six terms for G\ 3) in this way carries 
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(44) into 

ftWUu ft, -h, f P ) = bGHh)G l (-f 2 )G l (f p ) 

•[- 2 + tf u + + K u + - 2]. 
Going through all seven terms in Gi 2) carries (45) into 

^j<?i 2) (/i, /i, -/i, /„) = frGKWii-ftWiifpftKpiK* + 1) 

+ /?„(«„ + i) + (- 2) cur. + i) + ^(l) 

+ (-2) (-2) + (0)(tf u ) + (X.)(0)]. 

Substitution of the values of (?l°(/i, /2, — ft, f P ), I = 2, 3, 4, in the 
expression (42) for G 4 and combining terms leads to 

G<(fi, h, -h, U) = b'G\(h)G i {-h)G l Uv) 

-K U [2(K P + 1)(K,t 1) + Ki]. (47) 

When this is put in (40) we get the approximation we have been 
seeking for the third-order distortion (exp [j(2u>i — W2 + Up)Q) com- 
ponent of y(t). 

The procedure used to obtain (47) can also be used to show that 
the G4 in the expression (41) for the exp [j(a>i + a> p )£] component of 
y(t) has the value 

G*{f h f„ /» -U) = b*Gi(h)Q\Uv)Gi(-h) 

■Ku&iKt + 1)(K U + 1) + X p ]. (48) 

If we assume that the two series (40) and (41) converge at about the 
same rate, we can use (48) to get an idea of how large P p can be before 
the leading term in (■}()) ceases to be a good approximation to the 
typical third-order distortion term. Thus, we expect the leading term 
in (40) to be a good approximation as long as the ratio 

J|ppP,$Gi(/p)l 2 i:2(tfi + l)(Ku + 1) + K,3 (49) 

of the first two terms in (41) is small compared to unity. 

Note that setting ft = f\ and then interchanging /i and f p in the 
expression (47) for G 4 (/i, f\, — ft, f P ) carries it into the expression 
(48) for G«(/i, f P , f P , — f P ). This is to be expected since G 4 (fi, ft, f3, fi) 
is a symmetric function of /i, f%, fs, ft- 
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